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We clarify the role of crystalline electric field (CEF) in- 
duced magnetic anisotropy in the ground state and spin- 
wave spectrum of cubic rare-earth materials with dominating 
isotropic magnetic exchange interactions. In particular we 
study the hexaboride NdBg which is shown to exhibit strong 
spin-quadrupolar coupling. The CEF scheme is analyzed and 
a non-collinear magnetization response is found. The spin ori- 
entation in the antiferromagnetically ordered ground-state is 
identified. Moreover, the spin excitations are evaluated and in 
agreement with inelastic neutron scattering a suppression of 
one of the two magnetic modes in the strong-coupling regime 
is predicted. 

PACS numbers: 71.70.-d, 75.10.-b, 75.50.-y ,75.10.Jm, 
75.40.Gb, 78.70.Nx 



Introduction 

Over the last two decades cubic rare-earth hexaborides, 
REB6 (RE, rare-earth element), with CaBg-type crystal 
structure have been at the center of numerous studies 
of materials with crystallinc-clcctric-ficld (CEF) driven, 
non-trivial ordering phenomena. Among these com- 
pounds, CeB 6 (e.g., M) serves as a prototypical sys- 
tem which exhibits an impressively complex phase dia- 
gram. In this material the CEF of cubic symmetry selects 
the Tg quartet to be the ground state of the Ce 3+ ions 
[J = 5/2). The latter quartet is well separated from 
the next-highest T-? doublet by an energy gap of order 
540K |2|. Thus, on a low-energy scale, the physics of 
CeB 6 is reasonably well described by projecting onto the 
Tg subspace. Similar systems with Tg ground states can 
be realized starting from the right side of the rare-earth 
series, i.e., invoking compounds of cubic symmetry with 
Yb 3+ or Tm 2 ions, whose incomplete /-shell contains 13 
electron, or one /-hole. In accordance with Hund's rule 
and contrast to the Ce-case however, the Tg basis has to 
be constructed from a J = 7/2 multiplet, breaking direct 
electron-hole symmetry thereby. 

In this brief note we will focus on the hexaboride NdBg. 
Although investigated in detail experimentally by inelas- 
tic neutron scattering (INS) the anisotropy of the 
magnetically ordered state below the temperature Tc of 
order T c « 8.6K §,§ remains unclear as well as the exis- 
tence of only a single magnetic mode as observed by INS. 
The aim of our work is to consider these open issues. 



Crystalline electric field 

The CEF level scheme of the Nd 3+ multiplet (three 
/-electrons, J = 9/2, S = 3/2, L — 6) is consistent with 

the sequence T ( g\o K) - r^(135 K) - r 6 (278 K) [|. 
Similar to CeBg the energy gap separating the lowest 
quartet is large enough to restrict the Hilbert space to 
Tg ' only. The basis states of the latter sub-manifold can 
be represented by the linear combinations p|: 



^ +T = V! |+9/2) + v 2 |+l/2> + v 3 |-7/2) , 
=w x |+5/2) +w 2 |-3/2), 

tP +1 = Vl |-9/2) + v 2 |-l/2) + v 3 |+7/2) , 
=toi |-5/2) +w 2 |+3/2). 



(1) 
(2) 



The coefficients Vi and uii depend on the respective ratios 
of the CEF splittings. For NdB 6 one finds 



Vi = 0.1437, v 2 = -0.3615, v 3 = 0.9212; 



(3) 



wi = -0.9223, w 2 = 0.3865. 

The states in ([!]) and (||) have been labeled such that the 
2 nd index denotes a "spin"-like projections, whereas the 
1 st index stands for two "orbital" -like components which 
reflect the different shapes of the electron wave- functions. 
This leads to a description of the quartet in terms of two 
Pauli matrices, er and r p|||,[| 

0- z i> T ± = ±l/2Vv±, cr ± ^ T=F = 1p T ±, 
T Z lp±a = ±1/2V>±<7, T ± ^ a =^± cr . 

where the magnetic-moment operator can be represented 
in terms of er and r by: 

M a = [i B (£ + 2r]T a ) cr Q , (a = x,y,z). (4) 

Here T is a vector with components 

T x = ~2 Tz ~2 Tx ' ^ y ~ ~~2 Tz 2~ Tx ' ^ z = Tz ' (5) 

which transforms according to the T3 representation. 

Note, that for NdBg , £ and r\ depend on the parameters 
of CEF splitting @ with: 



i = -0.661 , 77 = -6.857 



(6) 



This identifies NdBg to be a system with strong cou- 
pling of the magnetic and quadrupolar degrees of free- 
dom. In the l/-electron(hole) cases £ and r\ are univer- 
sal and do not depend on the CEF splitting parameters. 
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For Ce 3+ £ = 2 and 77 = 8/7, for Yb 3 + £ = -8/3 and 
77 = —32/21. Therefore CeBg and YbBg exhibit rather 
weak spin-quadrupolar coupling with a characteristic pa- 
rameter 77/ (2f) = 2/7. 



Exchange anisotropy 

In this section we clarify the spin orientation in the 
magnetically ordered ground state. Most likely, the dom- 
inant interaction in NdBg is of isotropic magnetic ex- 
change type [Q. However, due to the T 8 ground state a 
CEF induced magnetic anisotropy exists which depends 
on the ratio £/rj. This can be understood easily by con- 
sidering the single-ion Zeeman interaction, i.e., — H • M 
in an external magnetic held H. The eigenvalues A in the 
Tg-space are 



A = 



±J^ +71 2 ± \ V \ { ^_ 2e) _ 3F{n){ ^_ 2e) , (7) 



measured in units of gfisH/2. This clearly manifests a 
cubic anisotropy through the function F(n): 



F(n) = n% 



n = H/H 



(8) 



The anisotropy results in a non-collinearity of the mag- 
netic field and the magnetization for any general orien- 
tation of H. Exceptions are the directions [111], [HO], 
[001], and their crystallographic equivalents. Energeti- 
cally favorable states are related either to the cubic axes 
([001]-type), if \q\ < 2|£|, or the cubic diagonals ([111]- 
type), if 1 77 1 > 2|£|. The anisotropy caused by the CEF 
disappears, if \r)\ = 2|£|. Therefore, we may conclude 
that Ce 3+ , Yb 2+ and Tm 3+ r 8 compounds tend to ex- 
hibit "easy axis" anisotropy (77/ (2^) = 2/7), whereas for 
Nd 3+ in NdBg we have 77/ (2£) rs 5.19 which results in 
"easy diagonal" anisotropy. 

Within a mean-field treatment of the exchange inter- 
action 



RR' a R 



'R' 



(9) 



R,R' 



where Jrr' is the exchange integral and Sr the spin at 
site R the magnetic field in (0) and (g) has to be replaced 
by the Weiss field Jo(S)/(<7/zb) with J = J2w 4r' if 
ferromagnetic exchange is dominant. The Lande factor 
in NdBg is g = 8/11. For bipartite antiferromagnetism 
(AFM) , the Weiss field on sublattice A is proportional to 
-J (Sa) + Ji(Sb) with J 0(1) = (-) £) R/ J RR - for R and 
R' on equal(opposite) sublattices. On sublattice B, one 
should replace A «-> B. 

Therefore, in conclusion, we expect [111] orientational 
ordering in the ground state of NdBg if isotropic exchange 
interactions are dominant IttOfl. 



Magnetic excitations 

In this section we focus on the spin dynamics by con- 
sidering the time-dependent magnetic susceptibility 



x s a p(k,t)=ie(t){[s ak (t),s^ k }) 



(10) 



Lower Greek indices of x an d the spin operator refer 
to x, y, z and boldfaced vectors k denote the momen- 
tum. We use a spin operator rescaled by 77 -1 , i.e., 
S a k = M a ic/(gfiBf]). Therefore the dependence of the 
magnetic spectrum on the CEF can be expressed solely 
in terms of the ratio £/r]. To evaluate ( [To|) we proceed via 
a mean field analysis consistent with AFM ordering || 
on a bipartite lattice. Rather than employing the Pauli- 
matrix representation @H§ of © we perform this anal- 
ysis using a dyadic basis to express the spin operator 
within the r§ manifold fTl"|| : 



N £ ' 



R 



R 



R' 



(11) 



where a summation over repeated indices is implied for 
the remainder of this paper and 



% = 



|A«R}<i/R| , b£, = |/jR')(z/R'| 



(12) 



are the dyades on sites R(R') of the magnetic A(B)- 
sublattice. \fi) are the eigenstates of the z-component 
of the spin in the Fg-manifold S a =z\fJ>) — s /j|m)- The 
spin should be quantized along (against) the [111] di- 
rection of the Weiss-field on the A(B)-sublattice sites. 
S£ v are the matrix elements of the spin correspond- 
ing to the latter quantization direction. The dyadic 
transition operators a£ and b£ v with /j, v = 1 ... 4 
can be recast into a 32-component operator A"^~ l "' 32 = 



{• 



(1,1),. ..(4,4) ,(1,1),.. .(4,4) 



} with a corresponding 32x32 



matrix-susceptibility of the AZ operators 



x ^(Kt) = ie(t)([A^t),A^]) 



(13) 



The original magnetic susceptibility (|o|) can be obtained 
from this by projecting the dyades onto the magnetic 
moment 



Xa/9(k,t) = i^"(k,t)C^ 



(14) 



(1,1),. ..(4,4) 



where Cf a = with ]r;^ = {S { <> 

fi ,(i,i),...(4,4)j j g 32-component vector for each spin com- 
ponent a. 

To proceed we evaluate the equation of motion (EQM) 
of the dyadic susceptibility 
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W"(M) 



-ie(t)([[AZ(t),H],A^}) 



(15) 



In this brief report we concentrate on the spin dynam- 
ics for next-neighbor (NN) AFM exchange-couplings J 
only. The effects of longer-ranged couplings will be dis- 
cussed elsewhere fl. Therefore, setting Jn 2 /g 2 to unity 
the Hamiltonian in terms of the dyades reads 



H = 



R,l 



O n J n Ur> U 



R "R+l 



(16) 



where 1 runs over the NN sites of R. The real-space rep- 
resentation of the commutator on the r.h.s of the EQM 
is evaluated using the algebra of the dyades, yielding 



i 



'R 



5ujfj, oXaj Xa 



(17) 



An analogous expression results on the B sub-lattice. On 
the mean-field level the EQMs are closed by factoriz- 
ing all quadratic terms in (|l7j) according to the scheme 



Or b w = («r>r' + Moreover, 'up'('down') 



[lll]-polarization on the A(B) sub-lattice is enforced by 
setting 



In momentum space the linearization results in 



(18) 



, H] = zS 4 a 4 (S^S^ - S^5 v °)ai° 



7 kl2 °k ) 



(19) 



where z is the coordination number and Z7k = Y]\ e • 
A similar equation arises for , H] introducing two ad- 
Switching to 



ditional 16x16 matrices L 



k22 



and L 



k21 



frequency space the EQMs can be solved as 
xf/3(k,w) =>1 



^Lk) 1 Xo c Po 



(20) 



where boldfaced symbols refer to matrix notation in a 
32x32 space. Lk is set by L^ Xa with i,j = 1,2 la- 
beling four 16x16 sub-blocks. Similarly Xo consists of 
four sub-blocks \o,i 3 wlth Xo,;#j = and Xo,n(22) = 

^^1(4)^1(4) _ ^^1(4)^1/1(4) _ 

Eqn. ( po|) allows for substantial simplifications. First, 
all diagonal dyades, i.e. a(b)^, commute with H. Sec- 
ond, the linearized form of ([17]) for the non diagonal 
dyades, i.e. for a(b)^ with /i ^= v, is diagonal with re- 
spect to /iv and remains local for nearly all pairs \xv. This 
follows from the identity pi 



S^S^ = 



(21) 



The only set of dyades which couple dispersively via the 
EQMs is 



(3,1) ,(3,4) ,(4,2), 
u k ' "k ' °k J 



(22) 



and the corresponding hermitian conjugate set _B^ _1 " 4 . 
From the preceding discussion it is conceivable that the 
complete spin dynamics can be expressed in terms of the 
physically relevant dyades B^T 1 "' 4 only. In fact, af- 
ter some elementary rearrangements of the matrix-EQM 
(pp|), the longitudinal spin susceptibility, which, due to 
cubic symmetry, is identical to the three-trace XaaO i J UJ ) 
simplifies to 



X s aa {^) = ~Tr[D- l N] 



(23) 



where the dynamical matrix D and the static susceptibi- 
lity-matrix N are identical to ((u>/z)l — Lk) and Xo^aa/z 
restricted to within the 4 dimensional subspace spanned 
by eqn. ([22]). The complex conjugate dyades B£ ^ intro- 
duce an overall prefactor of 2 only. After some algebra 
we find that D and N are determined by five parameters 
a, 6, c, d, and e through 



D = 



w — a — C7k — e7k 

w — b — e7k — G?7k 

C7k — e7k w + a 

—e7k G?7k w + b 



c — e c e 

e — d e d 

-c e —c —e 

e —d e d 



(24) 



with w = uj/z and 

a = S 44 (S'^ 1 — S^ 1 ) , b — S 44 (S 1 t 1 — S 133 ) 
c = >S a 5^ , d — S a S a , e — S a S a — V * cd . (25) 



J a a J 



With this the longitudinal spin susceptibility of eqn. (|2 
is obtained readily as 



xL( k >w) 



Z(k, w)/z 



(w 2 — wf)(w 2 — w 2 ) ' 
where the weight factor Z(k,w) given by 



(26) 



Z(k, w) = 2(ac -bd-(c- df^)w 2 

+ (ad - bc)(ab + (ad - bc)y k ) , (27) 

and the excitation energies iwi^k) are being set by the 
roots of the biquadratic equation 



w 4 + w 2 ((c-d) 2 lk -(a 2 + b 2 )) 

+a 2 b 2 - (ad - 6c) 2 7k = 



(28) 



In fig.(|l|) the dispersion as well as the weight i? 1-2 (k) = 
xf Q (k,a;)(w - Wi )2 (k)| it , =Wl2 ( k ) of the two positive-fre- 
quency modes is depicted along a path in the Brillouin 
zone (BZ) ranging from k = (1, 1, 1) to (0, 0, 0) to (1, 0, 0) 
for various values of the anisotropy ratio ^jr\. This figure 
clarifies the concluding issue aimed at in this brief note, 
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i.e. the observation of only a single excitation mode in 
NdB 6 . Based on the eigenvalues (0) two excitations of 
comparable energy are expected in the Weiss field of the 
AFM state at £/?y <C 1. However, fig. (Q) shows that 
only a single mode carries significant weight at small £/rj. 
We have shown this feature to remain valid for arbitrary 
range couplings Furthermore, in agreement with the 
spectrum of a single-ion pseudo-spin J — 3/2, the system 
exhibits a single-mode spin-wave like excitation at the 
isotropic point 2£ = rj. Only for intermediate anisotropy 
both modes show sizeable weight at any given point in 
the BZ. 



Conclusion 

In summary we have considered rare-earth compounds 
of cubic symmetry with a Ts-quartet ground-state of the 
RE ions. Particular emphasis has been put on the hexa- 
boride NdBg . Analyzing the CEF splitting we have iden- 
tified NdBg to be a genuine example of a system with 
strongly coupled magnetic and quadrupolar degrees of 
freedom. 

We have studied the CEF induced intrinsic magnetic 
anisotropy superimposed onto an isotropic exchange in- 
teraction revealing that NdBg should display magnetic 
anisotropy of a different type, i.e. 'easy diagonal', as 
compared to Ce or Yb compounds which show 'easy axis' 
anisotropy. 

The magnetic anisotropy leads to a non-collinear M vs. 
H behavior and is tempting to speculate that angular- 
dependent magnetization measurements on the corre- 
sponding RE cubic compounds, as well as diluted sys- 
tems, e.g, Lai-^Ce^Bg, should be able to detect this be- 
havior. 

We have evaluated the magnetic excitations in the 
AFM state of an 'easy diagonal' type using a dyadic oper- 
ator approach. For systems with strong spin-quadrupolar 
coupling this method is superior |j] to less controlled 
pseudo-particle descriptions which are applicable to the 
weak-coupling system CeB 6 @,||,U an d are based on the 
conventional <t-t Pauli-matrix representation (Q). In ac- 
cordance with the number of independent Pauli matrices 
(<r and r), we find two branches of spin excitations. How- 
ever, the spectral weights in the two magnetic channels 
are very different in a strongly coupled spin-quadrupolar 
system. In fact, in the £ = limit one channel disap- 
pears completely. This is reminiscent of the INS data on 
NdBg Q| which display only one branch of spin excita- 
tions. Although derived by a linearization of the EQMs 
we believe that our results are quite robust against non- 
linear corrections since the spin-wave spectrum in non- 
isotropic case is gapful. This should diminish the rele- 
vance of quantum fluctuations. 

Finally, regarding a direct comparison to experimental 
data we note that NdBg displays a [0, 0, 1/2] wave vector 
of the AFM modulation. This requires the inclusion of 




0.2 rt 



(1,1,1) k[7r] (0,0,0) (1,0,0) k[n] (0,0,0) (1,0,0) 
(1,1,1) 

FIG. 1. Dispersion and weight of spin excitations. 



longer-range exchange interactions which have been ne- 
glected in this paper. These and the physical nature of 
the branch unobservable by INS will be studied elsewhere 
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